#include<iostream>
#include "polynomial.h"
#include "vector"
#include <cmath>  
using namespace std;

int main(){
    vector <double> x(10), y(10);
    x[0] = x[1] = 0, x[2] = x[3] = 3, x[4] = x[5] = 5, x[6] = x[7] = 8, x[8] = x[9] = 13;
    y[0] = 0, y[2] = 225, y[4] = 383, y[6] = 623, y[8] = 993;
    y[1] = 75, y[3] = 77, y[5] = 80, y[7] = 74, y[9] = 72;
    Hermite_Interpolation p(x, y);
    Polynomial P = p.getPoly();
    cout << P << endl;
    cout << P.diff() << endl;
    cout << P(10) << ' ' << P.diff(10) << endl;
}